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Next-to-next-to-leading contributions to the cross section a{e^e~ T^r~) at energies close to threshold are 
analysed, taking into account the known non-relativistic effects and 0{a^) corrections. The numerical changes 
with respect to previous works are small, but the new corrections give a true estimate of the uncertainty in the 
theoretical calculation. 



1. Introduction 

The Tau-Charm Factory, a high-luminosity 
(~ lO'^^ cm~^ s~^) e+e" collider with a centre- 
of-mass energy near the t'^t~ production thresh- 
old, has been proposed as a powerful tool to per- 
form high-precision studies of the r lepton, charm 
hadrons and the charmonium system jl). In re- 
cent years, this energy region has been only par- 
tially explored by the Chinese BEBC machine 
(~ 10^"'^ cm~^ s^^). The possibility to operate the 
Cornell CESR collider around the T+r~ threshold 
|2j has revived again the interest on Tau-Charm 
Factory physics. Also the CMD-2 detector at the 
Novosibirsk VEPP-2M e+e^ collider will be soon 
ready to collect new data in this region 

A precise understanding of the e^e^ — > r+r^ 
production cross section near threshold is clearly 
required. The accurate experimental analysis of 
this observable could allow to improve the present 
measurement Q of the r lepton mass. The cross 
section a{e~^e~ —^ t~^t~) has already been anal- 
ysed to 0(a'^) in refs. including a resumma- 
tion of the leading Coulomb corrections. 

The recent development of non-relativistic ef- 
fective field theories of QED (NRQED) and QCD 
(NRQCD) H has allowed an extensive investiga- 
tion of the threshold production of heavy flavours 
at e'^e^ colliders. The threshold bb and tt 
js) production cross sections have been computed 
to the next-to-next-to-leading order (NNLO) in 
a combined expansion in powers of as and the 
fermion velocities. Making appropriate changes, 
those calculations can be easily applied to the 



study of T+r production One can then 

achieve a theoretical precision better than 0.1%. 

2. Perturbative calculation to ©(a*) 

At lowest order in QED, the r leptons are pro- 
duced by one-photon exchange in the s-channel, 
and the total cross section formula reads 

a.(e+e-^r+r-) = ^.(3-.2), (1) 

where v = -^/l — 4r7i^/s is the velocity of the final 
T leptons in the center-of-mass frame of the e+ e~ 
pair which makes as vanish when v —^ 0. 

Electromagnetic corrections of 0{a) arise from 
the interference between the tree level result and 
1-loop amphtudes. A factor a/v emerges in the 
1-loop final state interaction between the tau lep- 
tons, making the cross section at threshold finite. 
Furry's theorem guarantees that contributions to 
cr(e^e^ r+T^) coming from initial, intermedi- 
ate and final state corrections completely factor- 
ize at 0{a^), including real photon emission. 

Some undesirable features appear at 0{a'^): 
The two- loop T+T~7 vertex develops an c? jv^ 
term which makes the cross section ill-defined 
when w ^ 0, and multiple photon production of 
tau leptons by box-type diagrams and the non- 
zero interference of initial and final state radia- 
tion spoil exact factorization. However, as it has 
been recently shown [p^ , the squared amplitude 
of the e+e~ t^t^ box diagram (see Fig. Wj is 
proportional to a^v^, and so represents a N LO 
correction in the combined expansion in powers 



Figure 1. Direct and crossed box diagrams for 
e+e- ^ ff. 

of a and v, far beyond the scope of this analy- 
sis. In addition, contributions to the total cross 
section from diagrams with real photons emitted 
from the produced taus can be shown to begin at 
N'^LO, and factorization remains at NNLO. The 
total cross section can thus be written as an in- 
tegration over the product of separate pieces in- 
cluding initial, intermediate and final state cor- 
rections: 
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The radiation function F{s,w) ||Tl| describes ini- 
tial state radiation, including virtual corrections. 
The integration accounts for the effective energy 
loss due to photon emission from the e+ e~ pair, 
and it includes the largest corrections coming 
from the emission of an arbitrary number of ini- 
tial photons, which can sizable suppress the total 
cross section, a collects only final-state interac- 
tions between the tau leptons, and it is usually 
written in terms of the tau spectral density R^, 
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with a-pt — . The threshold behaviour of the 
total cross section will be ruled by the expansion 
of at low velocities. 

3. Non-Relativistic Corrections: NRQED 

The aim of the effective field theory approach 
to threshold particle production is to achieve a 
given accuracy in the combined expansion in pow- 
ers of a and the velocity v. Such double expansion 
is needed at threshold due to the appearance of 
(a/w)" terms in the QED perturbative expansion 
at the n-th loop order, which forces one to treat 
a and v on the same footing. A non-perturbative 



Figure 2. Ladder exchange of photons between 
produced fermions. 

procedure to deal with such singular terms in the 
limit w ^ is therefore mandatory. 

The leading divergences (i.e. (^) , n > 1), 
emerging from the ladder diagrams with Coulomb 
photons shown in Fig. ^, are resummed in the 
well-known Sommerfeld factor |l2| 
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multiplying the Born cross section Next-to- 
leading order (NLO) entails terms proportional 
to (a/u)" X [a,v], while NNLO accuracy stands 
for contributions (a/w)" x [a'^,v'^,av]. 

A systematic way to calculate these higher- 
order corrections in this regime requires the use 
of a simplified theory which keeps the relevant 
physics at the scale Mv ~ Ma, characteristic of 
the Coulomb interaction. NRQED |^ was de- 
signed precisely for this purpose. It is an effective 
field theory of QED at low energies, applicable 
to fermions in non-relativistic regimes, i.e. with 
typical momenta p/M ^ v <C 1. Interactions 
contained in the NRQED Lagrangian have a def- 
inite velocity counting but propagators and loop 
integrations can also generate powers of v. With 
appropriate counting rules at hand, one can prove 
that all interactions between the non-relativistic 
pair T+T~ can be described up to NNLO in terms 
of time-independent potentials jl^ , derived from 
the low-energy Lagrangian. Therefore, the low- 
energy expression of the r spectral density is re- 
lated with the non-relativistic Green's functions 

= ^ GiE) - G^iE)) (5) 

with Ci a short distance coefficient to be deter- 
mined by matching full and effective theory re- 



suits and E = rrirV^ the non-relativistic energy. 
The details of this derivation can be found in the 
Appendix B of ^ . The Green's function G obeys 
the Schrodinger equation corresponding to a two- 
body system interacting through potentials de- 
rived from £nrqed at NNLO: 

X G{r,r',E) = S^^\r ^ r') (6) 

The term — 4^p- is the first relativistic correction 
to the kinetic energy. Vc stands for the Coulomb 
potential with ©(a^) corrections. At NNLO, the 
heavy leptons are only produced in triplet S-wave 
states, so we just need to consider the corre- 
sponding projection of the Breit-Fermi potential 
Vbf- Finally, Va^ is a NNLO piece derived from 
a contact term in £nrqed, which reproduces the 
QED tree level s-channel diagram for the process 

A solution of eq. (|^) must rely on numerical 
or perturbative techniques. In the QED case, 
a significant difference between both approaches 
is not expected, being a such a small parame- 
ter. Consequently wc follow the perturbative ap- 
proach, using recent results where NLO and 
NNLO corrections to the Green's function arc 
calculated analytically, via Rayleigh-Schrodinger 
time- independent perturbation theory around the 
known LO Coulomb Green's function Gc- We re- 
fer the reader to Appendix C of [D for complete 
expressions of the Green's function corrections. 

4. Vacuum Polarization 

For a complete NNLO description of a{e~^e~ 
T+T^), two-loop corrections to the photon prop- 
agator should be included. The light lepton con- 
tributions to the vacuum polarization are the 
standard 1- and 2-loop perturbative expressions. 
For the r contribution in the threshold vicinity 

> 4M^, resummation of singular terms in the 
limit V ^ is again mandatory. Under the as- 
sumption a ~ u, we need to know NLO contribu- 
tions to Ilr{q^), performing the direct matching 
for both real and imaginary parts. 

In the hadronic sector, we can relate the 
hadronic vacuum polarization with the total cross 



section a{e^e~ ^ 7* ^ had). Below 1 GeV, the 
electromagnetic production of hadrons is domi- 
nated by the p resonance and its decay to two 
charged pions. The photon mediated tt+tt" pro- 
duction cross section is driven by the pion elec- 
tromagnetic form factor F{s). An analytic ex- 
pression for F{s) was obtained in Ref. using 
Resonance Chiral Theory and the restrictions im- 
posed by analyticity and unitarity. The obtained 
F{s) provides an excellent description of exper- 
imental data up to energies of the order Sp ~ 
1 GeV2. 

For the integration region above Sp, we use 
Imllh^j as calculated from pQCD. Our simple es- 
timate has been proved ^ to deviate by less than 
5% for the running of a at the scale -y/s = Mz- 
Considering that flh^d modifies a{e'^e^ — > r+r^) 
near threshold by roughly 1%, our result has a 
global uncertainty smaller than 0.1% for the total 
cross section. Clearly, our antsaz could be easily 
improved using a more realistic hadronic spec- 
trum, and, indeed, most experiments have their 
own routines to accurately implement vacuum po- 
larization together with initial state radiation in 
their data analyses. 

5. Electroweak and Bound States efTect 

Electroweak production of heavy quarks in- 
cluding threshold effects has already been stud- 
ied in previous papers p^ , and can be easily in- 
corporated in our basic formula (^. However, 
the characteristic m^/A/f suppression of this set 
of corrections and the small value of the neutral- 
current couplings of the leptons make these terms 
fully negligible in our analysis. 

Green's functions develop energy poles below 
threshold corresponding to spin triplet (n^S*!) 
electromagnetic r+r" bound states. The small 
width of these bound levels, dominated by their 
e+e~ decay rate 

_ rrira^ _ 6.1 • 10"^ 

make these states very difficult to be resolved ex- 
perimentally. For the same reason, the bound 
states cannot affect the shape of the cross section 
at threshold, contrary to the case of heavy quark 
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Figure 3. The spectral density Rt at low veloci- 
ties in both QED and NRQED. 

threshold production, where bound states play a 
crucial role. 

6. Numerical analysis for (T(e+e^ t^t^) 

The need for performing resumma- 
tions of the leading non-relativistic terms 
{a/v)" [v,va,v^, . . .] is evidenced in Fig. ^. The 
QED spectral density vanishes as u ^ 0, due to 
the phase space velocity in formula (|l]), which 
is canceled by the first term appearing in 
the 0{a) correction, making the cross section 
at threshold finite. More singular terms near 
threshold, . . . arising in higher-order correc- 
tions completely spoil the expected good conver- 
gence of the QED perturbative series at low v. 
This is no longer the case for the effective theory 
perturbative series, whose convergence improves 
as we approach the threshold point, as shown in 
Fig. ^, and higher-order corrections reduce the 
perturbative uncertainty. In the whole energy 
range displayed in Fig. ^, the differences between 
the NNLO, NLO and LO results are below 0.8%, 
which indicates that the LO result, i.e. the Som- 
merfeld factor, contains the relevant physics to 
describe the threshold region. Adding the inter- 
mediate and initial state corrections we have a 
complete description of the total cross section of 
T~^T~ production, as shown in Fig. ^. Coulomb 
interaction between the produced r's, becomes 
essential within few MeV's above the threshold, 
and its effects have to be taken into account to 
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Figure 4. Relative sizes of corrections to Rr{s) as 
calculated in NRQED. 

all orders. Initial state radiation effectively re- 
duces the available center-of-mass energy for r 
production, lowering in this way the total cross 
section. 

We should emphasize that NNLO corrections 
do not modify the predicted behaviour of the 
LO and NLO cross section, but are essential 
to guarantee that the truncated perturbative se- 
ries at NLO gets small corrections from higher- 
order terms. Differences with previous evalu- 
ations mainly concern the correct resummation 
of NNLO terms ex (a/v)" x [a^,v^,av], which 
on perturbative grounds means that our results 
should differ by terms. The latter are seem- 
ingly negligible, but one has to check that all- 
order resummation of them does not produce any 
unexpected enhancement. 

In addition, factorized formulas proposed in 
previous works do not have under control the 
NNLO terms being resummed beyond the order 
in a at which the matching between the non- 
relativistic and the QED results is performed. 
These fake resummations may be harmless in 
QED, being a such a small parameter, but must 
be handled with care when estimating the un- 
certainties of the theoretical calculation. In that 
sense, notice in Fig. ^ that NNLO corrections 
to the spectral density are not 0{a^) ^ 0.005% 
with respect to the leading-order, as we would 
naively estimate. Hence, we can safely conclude 
that the theoretical uncertainty of our analysis of 
cr(e+e~ T+T^) at energies close to threshold 
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Figure 5. Total cross section (T(e+e^ t^t^) 
at threshold. 

is lower than 0.1%. 

A FORTRAN code which evaluates the spec- 
tral density i?™'"°(s) will be soon available at 
http:/ /alpha. ific.uv.es/^ruiz 
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